Itinerant ferromagnetism in half- metallic C0S2 
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We have investigated electronic and magnetic properties of the pyrite-type C0S2 using the linearized 
mufhn-tin orbital (LMTO) band method. We have obtained the ferromagnetic ground state with 
nearly half-metallic nature. The half-metallic stability is studied by using the fixed spin moment 
method. The non-negligible orbital magnetic moment of Co 3d electrons is obtained as fiL = 0.06/iB 
in the local spin density approximation (LSD A). The calculated ratio of the orbital to spin angular 
momenta {Lz)/{Sz) — 0.15 is consistent with experiment. The effect of the Coulomb correlation 
between Co 3d electrons is also explored with the LSDA + U method. The Coulomb correlation at 
Co sites is not so large, C/ < 1 eV, and so C0S2 is possibly categorized as an itinerant ferromagnet. 
It is found that the observed electronic and magnetic behaviors of C0S2 can be described better by 
the LSDA than by the LSDA -I- U. 

PACS number: 75.10.Lp, 71.15.Mb, 71.20.Be, 71.70.Ej 



I Transition metal disulfides of MS2 (M = Fe, Co, Ni, 
etc.) with the pyrite structure shows a variety of phys- 
ical properties. Electronic and magnetic properties of 
the MS2 has been qualitatively explained by the succes- 
sive filling of the Cg band of metal ions. In the NaCl- 

I type structure of pyrite MS2 , divalent metal ions of 

are located at the center of octahedron formed by 82^^ 

' dimers. Hence, the M-3c? bands are split into t2g and Cg 
subbands due to large crystal field. The M-3d electron 

■ configuration varies from ^23^"' ^%^g' ^Zg^g ^'^^ FeS2, 
C0S2, and NiS2. All of them are in the low-spin states, 
S- = 0, 5* = 1/2, and 5* = 1 for FeS2, C0S2, and NiS2, re- 
spectively. With fully occupied t2g and completely empty 
bands, FeS2 is a semiconductor with an energy gap of 
Eg ~ 0.8 eV M. C0S2 is a ferromagnetic metal with 
Tc — 120 K mM- NiS2 is an antiferromagnetic insulator 

I (Tjv — 40 K) with a half-filled band due to the large 
on-site Coulomb correlation 

] Despite extensive studies on these systems, the elec- 
tronic and magnetic properties of C0S2 are not fully un- 
derstood yet. The issues for C0S2 to be resolved are {i) 
whether it is a strongly correlated electron system, (m) 
whether it is classified as an itinerant electron ferromag- 
net, and (Hi) whether it is a half-metal or not. It is con- 
sidered that Co 3d states in C0S2 are just at the bound- 

I ary of the localized and delocalized regime |^ . Magnetic 
measurement using the polarized neutron diffraction in- 
dicates that most part of the magnetic moment is lo- 
calized at Co sites Q. The analysis of photoemission 
spectroscopy (PES) data for C0S2 in terms of the cluster 

, model calculations yields a large Coulomb correlation pa- 
rameter, U — 3.0 ^ 4.2 eV [^,0. According to the mag- 
netic circular dichroism (MCD) measurements, a non- 
negligible orbital magnetic moment is observed for Co 3d 
electrons, {L,)/{S^) = 0.18 § and 0.14 §. Such an un- 
quenched orbital moment is exceptional, considering the 
low-spin state of Co 3d electrons in C0S2 with a large 
crystal field splitting of lODg ~ 2.5 eV. On the other 
hand, transport properties of C0S2 are known to be well 
described by the itinerant electron model ||9|,p^. Further, 



the dynamical susceptibility measured by the inelastic 
magnetic neutron scattering pT[ | and optical properties 
|]l2| are analyzed based on the itinerant band model. 

Using the linearized atomic orbitals (LCAO) band 
method, Zhao et al. [|l3| has found that C0S2 is close to 
half-metallic in the sense that the occupation number of 
the minority spin Cg band is very small compared to that 
of the majority spin Cg band. In fact, recent experimen- 
tal reflectivity data have been consistently interpreted 
under the assumption of half- metallic C0S2 . In con- 
trast, the linearized mufhn-tin orbital (LMTO) band cal- 
culation by Yamada et al. gives rise to the normal- 
metallic ferromagnetic ground state with a partially filled 
minority spin Cg band. 

To clarify the above mentioned issues and to resolve 
the differences between existing band calculations, we 
have reexamined electronic and magnetic properties of 
C0S2 using the LMTO band method. Both the local spin 
density approximation (LSDA) and the LSDA -I- U ap- 
proximation |l^ are employed to explore the effect of the 
Coulomb correlation interaction between Co 3d electrons. 
The von Barth-Hedin form of the exchange-correlation 
potential is utilized. Further, the fixed spin moment 
method is used to check the stability of half-metallic 
state. We have also estimated the orbital contribution 
to the total magnetic moment in C0S2, using the fully 
relativistic calculations in which the spin-orbit coupling 
is simultaneously included in the self-consistent varia- 
tional loop |l^. Two main parameters in the LSDA -|- 
U method are the Coulomb U and the exchange J inter- 
actions. The relations of these parameters to the Slater 
integrals are given hyU = and J = {F'^+F^)/1A. The 
ratio of F'^/F'^ is known to be constant around 0.625 for 
most 3d transition metal atoms [0. We have used pa- 
rameter values oi U ~ 0.0 ^ 4.0 eV with a fixed J — 0.89 
eV. 

The crystal structure of pyrite C0S2 is cubic with the 
lattice constant of a = 5.528 A and the space group 
is T^{Pa3). There are twenty-four symmetry operations. 
The four Co atoms are located at positions (4a) : (0,0,0), 
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FIG. 1. The LSDA DOS of pyrite C0S2 at each atomic 
site. The two peaks at the highest binding energy side (upper 
panel) are bonding and antibonding S 3s bands, respectively. 
The antibonding S Spa* band is located at about 2.5 eV above 
Ep. The spin split Co 3d bands, t2g and Eg, axe formed be- 
tween S 3p and S Spa* bands (lower panel). The eg band near 
Ef manifests nearly half-metallic nature of C0S2. 

(1/2,0,1/2), (0,1/2,1/2), and (1/2,1/2,0). The eight S 
atoms are at positions (8c) : {u,u,u), (u4-l/2, u, w-f 1/2), 
(m, u-I-1/2, M-f 1/2), and (w-f 1/2, u+1/2, u) where we have 
used the position parameter of m = ±0.389. To improve 
the packing ratio, twenty-four empty spheres are intro- 
duced at positions {24d). For the Brillouin zone (BZ) in- 
tegration, the tetrahedron method is adopted with 80 k- 
points sampling in the irreducible BZ wedge. The LMTO 
basis functions are included up to / = 2 for Co, and I — 1 
for S and empty spheres. 

In Fig. |l|, we show the LSDA density of states (DOS) 
obtained at the experimental lattice constant. Because 
of strong covalent bonding nature of S-S dimers, each S 
band can be identified with molecular orbitals. The two 
peaks in the highest binding energy side are the bonding 
S 3sa and the antibonding S 3sa* bands, respectively 
. The broad band of about 6 eV width between —4 
and —10 eV consists mainly of S 3p, which is a mixture 
of Spa, 3p7r, and Spvr* bands. The antibonding S Spa* 
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FIG. 2. The LSDA total energy as a function of lattice 
constant. The arrow indicates the experimental lattice con- 
stant. The ferromagnetic ground state is lower in energy than 
the paramagnetic state far below the experimental lattice con- 
stant. 

intra-dimer band is located at about 2.5 eV above the 
Fermi energy Ep. The spin split bands of Co t2g are 
located between —1 and —3 eV below Ep while Co Cg 
band are near Ep. It is seen that Co 3d bands are formed 
in the energy range between S 3p and S Spa* bands. 
The t2g band is nearly dispersionless with very narrow 
band width, whereas the eg band near Ep is dispersive 
with relatively large band width of about 2.5 eV due 
to strong hybridization with S 3pa* band. Nearly half- 
metallic nature and the exchange splitting Ax ~ 1.0 eV 
are prominent for Cg band. 

The present results are slightly different from those of 
existing calculations (l^Jl^ • In Rcf. S Spa* band is at 

1.5 eV above Ep which is lower than the present result 
by about 1.0 eV. Whereas, in Ref. the states of S 3pa* 
band is smeared out by too strong hybridization with 
Co Cg bands. The bremsstrahlung isochromat spectrum 
(BIS) shows clear double peak structure at about 1.0 and 
2.5 eV above E^ which is consistent better with 

the present result. That is, the 1.0 eV peak corresponds 
to the minority spin eg band and the 2.5 eV peak to 
S Spa* band. The occupation of the minority spin Cg 
band is only 0.04 electrons/Co, reflecting that C0S2 is 
almost half- metallic. This feature is in agreement with 
that obtained by Zhao et al. In the analysis of the 
reflectivity data in terms of the half-metallic C0S2 , 
the minority spin Cg band is assigned at about 1.5 eV 
above Ep. This is close to but a bit higher than the 
present LSDA result. 

Figure ^ is the LSDA total energy as a function of the 
lattice constant. The ferromagnetic ground state is lower 
in energy than the paramagnetic state even far below 
the experimental lattice constant of a = 5.528 A repre- 
sented by the arrow. The minimum of the total energy 
is located near the experimental lattice constant for both 
the paramagnetic and the ferromagnetic states, and thus 
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FIXED SPIN MOMENT (Hb) 

FIG. 3. The behavior of the spin magnetic moment at 
each site (upper panel) and the total energy (lower panel) 
with varying the total spin magnetic moment M{^b) in the 
unit cell. The arrow denotes the normal-metallic ferromag- 
netic ground state with M = 3.66/ifl. The half-metallic state 
with integer M = is very close in energy to the nor- 

mal-metallic ground state. 



the LSDA seems to describes well the cohesive bonding 
properties of C0S2. The ferromagnetic ground states is 
lower in energy by about 10 mRy than the paramagnetic 
state at the experimental lattice constant. The energy 
difference between the ferromagnetic and the paramag- 
netic states decreases with decreasing the lattice con- 
stant. This is contrary to the result of Yamada et al. 
in which the paramagnetic is stable near the experimen- 
tal lattice constant and the ferromagnetic state becomes 
stable only for a > 5.57 A. As they pointed out, it seems 
to be ascribed to small number of k-points they used. 
On the basis of the total energy curve shown in Fig. ||, 
the observed metamagnetic transition behavior under the 
high magnetic filed, the high external pressure, and the 
chemical pressure by substituting Se in place of S can be 
qualitatively understood |10 2^ 



TABLE I. Calculated spin and orbital magnetic moments, 
/is and /iL (/is) of C0S2 for various Coulomb correlation U 
(eV). All the results are from the relativistic calculations. 
MfaifJ^s) is the total magnetic moment per formula unit. The 
total magnetic moment per Co atom utot = 0.86/ib and the 
ratio of {Lz)/{Sz) = 0.15 in the LSDA is consistent with ex- 
perimental data. 
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By using the fixed spin moment method within the 
LSDA, the total energy is evaluated with varying the 
magnetic moment in C0S2 (Fig. ||). For all the fixed spin 
moment, O.OO/is < M < 4.00/iB, the ferromagnetic state 



is stable as compared to the paramagnetic state. The 
spin direction of S atoms is always parallel to that of Co. 
The ground state is found for M = 3.66/is (denoted by 
the arrow) which is a normal ferromagnetic metal and 
lower in energy by about 10 mRy than the paramag- 
netic state. It is noticeable that the half-metallic state 
with exact integer magnetic moment of M = 4.00/iB is 
very close in energy to the normal-metallic ground state. 
The total energy difference between the normal-metallic 
ground state and the half- metallic state is less than 0.1 
mRy because of quite small electron occupation of the 
minority spin Cg band in the ground state. Thus it is 
expected that, with only small number of hole doping, 
the half- metallic state can be realized in C0S2 fl^ . 

In general, a polycrystalline half-metallic system is ex- 
pected to show large magnetoresistance (MR) behavior. 
It arises from the activated electron tunneling mecha- 
nism at the grain boundary with the applied external 
magnetic field ||2^. To our knowledge, there has been 
no report available on the MR measurement for C0S2. 
Therefore, the MR measurement on polycrystalline C0S2 
and Fea:Coi_a;S2 is desirable to exploit the correlation be- 
tween the half-metallic nature and the MR. 

We have applied the LSDA and the LSDA-I-J7 meth- 
ods to explore the Coulomb correlation effect on the 
magnetic properties. The observed magnetic moments 
per Co atom are in the range of 0.84 ^ 0.91/is [^,^. 
The orbital magnetic moment in C0S2 is obtained by 
including the spin-orbit coupling in the total Hamilto- 
nian. In the LSDA -I- U calculations, we have used 
fixed J(= 0.89 eV) which is insensitive to local environ- 
ment. In the LSDA, we have obtained the spin and or- 
bital magnetic moments per Co atom, /ig = 0.80/iB and 
liL = 0.06/iB, respectively (Table ^. This produces the 
ratio of {Lz)/{Sz) = 0.15, which is close to the MCD 
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measurement of {L^)/{S^) = 0.18 § and 0.14 g. The 
total magnetic moment per Co atom utot = 0.86/iB in 
the LSDA is also in agreement with the experiment. 

Both the spin and orbital magnetic moments increase 
with increasing the Coulomb correlation parameter U. 
For U }t 3.0 eV, the minority spin Cg band moves to 
higher energy states and C0S2 becomes completely half- 
metallic. For U = 4.0 eV, the minority spin eg band is 
formed at about 2.0 eV above Ep. However, the discrep- 
ancy of the total magnetic moment between the experi- 
ment and the theory becomes significant for U > 2.0 eV. 
Moreover, the ratio of {L^)/{S^) = 0.31 for U = 2.0 eV is 
too large compared to the experimental values ||^,^ . This 
finding suggests that the size of the on-site Coulomb cor- 
relation in C0S2 is to be rather small U "^1 eV, and the 
LSDA gives a better description of Co 'id states in C0S2. 
Note that the LDA results for FeS2 show an excellent 
agreement with PES and BIS data [^,|^. In addition, 
the weak satellite feature in the PES of C0S2 (|] seems 
to be consistent with the present results of rather small 
Coulomb correlation U . 

Finally, let us remark on the electronic structure of 
C0S2 in its paramagnetic phase (T > Tq). As mentioned 
in the introduction, it is an unresolved issue. The ques- 
tion is whether the exchange splitting in Co 3c? bands 
persists above Tc or not. The present study reveals that 
the width of Co Cg band is larger than both the Coulomb 
correlation interaction U and the exchange splitting Ax . 
It thus suggests that the Co 3d electrons near E-p will be- 
have as itinerant to manifest Stoner-type magnetic phe- 
nomena at the finite temperature. Indeed, there are such 
evidences from several experiments. Resonance PES data 
for Co 3c? bands taken at temperatures across the mag- 
netic transition indicate a slight but noticeable spectral 
change between the paramagnetic and the ferromagnetic 
phases, reflecting the long range spin exchange effect . 
Optical spectra above Tq are described well using the 
itinerant Stoner model p2[ . Also inelastic magnetic neu- 
tron scattering data above Tc exhibit the Stoner excita- 
tions stemming from the itinerant electrons ||ll|. On the 
other hand, based on the PES studies, it is claimed that 
the electronic structure in the paramagne tic p hase should 
be described by the local band picture |||,|l^. Hence fur- 
ther studies are necessary to ascertain this point. 

In conclusion, we have investigated the electronic and 
magnetic properties of the pyrite-type C0S2 using the 
LMTO band method. We have found the following re- 
sults: (j) C0S2 does not belong to the strongly correlated 
electron system, because the effect of the Coulomb corre- 
lation at Co sites is rather small, U ^1 eV, (m) C0S2 can 
be categorized as an itinerant electron ferromagnet, since 
the band width of Co Cg state is larger than U and the 
exchange splitting Ax , {in) C0S2 is nearly half- metallic 
in that the half-metallic state is very close in energy to 
the normal-metallic ground state. A non-negligible or- 
bital magnetic moment of Co 3c? electrons is obtained 
as = 0.06^B in the LSDA, which is consistent with 
experiment. Accordingly, the electronic and magnetic 



properties of C0S2 are described better with the LSDA 
than with the LSDA-hC/. 
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